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Abstract. We study the dynamics of a spin coupled to an oscillating magnetic field, 
in the presence of decoherence and dissipation. In this context we solve the master 
equation for the Landau-Zener problem, both in the unitary and in the irreversible 
case. We show that a single spin can be "magnetized" in the direction parallel to the 
oscillating bias. When decay from upper to lower level is taken into account, hysteretic 
behavior is obtained. 
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1. Introduction 

This paper concerns the dynamics of a driven two-level system with decoherence and 
dissipation. Such a simple quantum-mechanical problem has been studied in various 
contexts, as, e.g., nuclear and atomic physics. What follows is predominantly motivated 
by molecular magnets and quantum information, as long as these topics suggest new 
perspectives to the old problem. 

The perspective of a comprehensive theory of molecular magnets stimulated a large 
amount of work. Rather than giving an extensive review, we merely quote few topics in 
this area. A 1/2 spin (the central spin) has been introduced in modeling the macroscopic 
spin tunneling (for a short review see, e.g., Ref. [|l|); dissipative spin reversal, incoherent 
Zener tunneling, decoherence were also considered |^ ^, |14|, ||, ||. 

Work more directly related to ours can be found in Refs. |0, ||, H- In classifying 
the behaviors of our system, we find effects such as single spin dynamical magnetization 
and hysteresis; the meaning of these terms will be made clear in the sequel. 

In the presence of an oscillating external field, the evolution is characterized by 
a sequence of level-crossings, where the Landau-Zener (L-Z) model |T0[ applies, 
separated by intervals of "normal" evolution. 

We classify the dynamics from the behavior in the L-Z regime and in the "normal" 
regime. In order to do that, we first solve the Landau-Zener problem in the master 
equation for the density matrix. This is done both in the unitary and in the irreversible 
case. 

Our starting point is Kraus's representation theorem, for the density matrix of 
a two-level system: a short reminder of the theory is given in Sections 2,3; for the sake 
of completeness, in Appendix I we recall the conditions for the complete positivity of 
the evolution ||T2| . 

In Section 4 we discuss the unitary case by means of qualitative methods. We 
examine high frequency fields as well as quasistatic fields. We show that a "magnetized" 
spin can be sustained, in the direction of a zero-average oscillating bias. 

The non-unitary case, where dissipation and decoherence are taken into account, is 
discussed in Section 5. We show that the L-Z problem can be solved in the irreversible 
case; the solution is explicit for two-level systems. 

We find the lowest order (in the 1/t expansion) renormalizations of the dissipation 
and decoherence parameters: the non-unitary processes do not modify the scattering 
properties of the solution at the level crossings. After analyzing the single-spin 
"magnetization" with decoherence, we take into account the relaxation from the upper 
to the lower level. In section 6 we show that an oscillating external field gives rise to 
hysteretic behavior. 
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2. Master equation for the two-level system 



We start with the following master equation, for the density matrix of a two-level system: 

3 
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-t[H, p{t)] + i E Pit) a]] + [a,p(t), a]]) 



(1) 



where H is associated with reversible dynamics, the matrix A is Hermitean and the CTj's 
are the Pauli matrices. 

The Hamiltonian is parametrized by 3 real parameters: 



and further 9 real parameters label the 3x3 hermitean matrix Aij. 



(2) 



Notice that this is consistent with the Kraus representation |13|, where the time 
evolution of a two-level system is characterized by 12 real parameters. 

As it is well known, the coefficients of A are the correlation functions of the 
macrosystem coupled with the spin. 

We then write Eq. ^ in the standard representation: 
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we have: 

v(t) = Mv(t) + C 
The unitary contribution in M has the form: 
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and the decoherence and dissipation processes give: 
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In Eq. H, the explicit form of M is made more readable by letting: 



71 = 2 (A22 + A33) 72 = 2 (A33 + An) 73 = 2 (An + A22) 
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3. Geometrical interpretation of the 12 parameters 

The phase space of the system is a Poincare surface, undergoing elementary 
instantaneous transformations, according with M and C. 

More specifically, while the antisymmetric part determines rotations, the symmetric 
part gives dilatations along 3 principal axes. The inhomogeneous term translates the 
surface. 

The connection with Kraus's evolution of the density matrix over a finite time 
interval is readily verified in the case of constant coefficients; in such a case, by explicit 
time integration, one gets the affine [0, [l^, map: 



(7) 



which again depends on 12 parameters. 

In our problem problem the hamiltonian is time-dependent: the Poincare surface 
still undergoes an affine transformation under a finite time evolution, but the operators 
Ai are fairly complicate functions of M[t) and Ci{t). 

We start by assuming a time-independent matrix A in eq. |l|. The eigenvectors of 
D identify the dilatation axes of the Poincare surface. The choice of representation for 
the hamiltonian part, given by M, fixes a second, time-dependent frame. 

Everything greatly simplifies if the two frames coincide, i.e. if D is diagonal: 



a 
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The parameters are then reduced from 12 to 9: 3 of them come from the 
hamiltonian, 3 from decoherence and dissipation, (71, 72, 73) , and the last 3 from the 
inhomogeneous terms Cj. 

E.g., in the case of an isotropic thermal bath with average phonon number n one 



has 16 



7i 



72 



In 



Ci — C2 



73 = -7(2n + 1) 
Cs = -7 



The Landau-Zener Model with Decoherence: 



The case S=l/2 



5 



4. Unitary evolution 

In this section we illustrate the dynamics in the unitary case. The external field is 
taken, without loss of generality, in the x-z plane: this implies A' = 0. We consider an 
oscillating z-component: = uJo{t), and a constant x-component: = A. 

The evolution can be depicted as a sequence of L-Z crossings and "normal" regimes. 

We characterize the "normal" regime with the condition |co'o(^)| ^ it is 
dominated by rotations in the plane X, Y, as one can verify by inspection of eq. |^. 

One must not expect that the representative on the Poincare surface retraces its 
path backwards in coincidence with the external field: in fact the solution is the time- 
ordered exponential of an operator. This is particularly true when co'o(t) is of the order 
of A. On the other hand, as long as the magnetic field {A, co'o(^)} runs over a line in 
the plane {B^, B^}, no Berry phase is to be expected. 

We take, as an initial condition, the "fully magnetized" configuration Z = ±1. We 
have then: 

X = -uJoY Y = uJoX - AZ Z = AY (8) 



ujo = Bocos{Qot) 

where Qq is the bias frequency. 

The Schrodinger equation corresponding to system |] has the form: 

idi = ^(tuofli + Aas) (9) 
ia2 = -(-LjQa2 + Aai) 

where Oj, (i = 1, 2) are the spinor amplitudes. 

In the atomic interpretation of the two-level system, A is the interlevel transition 
amplitude and uq the level spacing. 

In the Landau-Zener (L-Z) regime, i.e. close to the level crossing, assumed at t = 0, 
one has: ujq = Qt = B^Qot. The asymptotics of the solution, when ujq is linear in t, is a 
textbook topic; it can be readily obtained by WKB methods. Let us write the equation 
for the amplitude 02 (t), in the form: 

^Mt)= -W\t).a2{t) (10) 
W{t) = 

From the large time behavior of equation |10|, we get the semiclassical momentum 

^ ^ Qt u A2 
pit) = h - u = — 
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It is natural to choose the following scattering states: ea;p{±if2t^/4} ■ {t)^'^". 
The scattering matrix 5* is obtained from the Weber functions, who solve Eq. |10| (see, 



e.g., [T7|, and references therein). The procedure is straightforward once one realizes 



that eq. [ly is the Schrodinger equation for an inverted harmonic oscillator m. D = 1. 

The 5* matrix identifies a rotation. The branching properties of the solution, 
characterized by the single parameter z/, uniquely determine S and the rotation. Since 
z/ is a constant in our system, each crossing produces an identical rotation. 

Let us label S with 0, ch: 



2tt 
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= cos{9) = 
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= i sin{0) ■ 


exp{—i(j)} 
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where r(x) is the Gamma function. Clearly Q acts on the population Z, and on the 
coherences X, Y . 

One finds, for the ratio of "in" and "out" populations, (T = Z(+oo)/Z(— 00)): 
T = T{y) = 2 ■ exp{-27iu} - 1. 

So far, we discussed the Schrodinger equation. Here we solve the corresponding 
master equation problem. Notice that, since in going from the spin variables ai,2 to 
the vector variables Vi of the density matrix: Vi = a*ia\'^ak {i = 1,2,3), {l,k = 1,2), 
one performs a nonlinear transformation, the new problem is not a straightforward 
translation of the old one. 

Rather than satisfying a second-order order equation as in Eq. |10|, each vector 
component satisfies a third-order equation, so that the connection with the Schrodinger 
equation is lost. 

In spite of that, in the representation of frequencies u, the population Z^u) satisfies 
a Schrodinger equation for an (inverted) harmonic oscillator in D=2, with a centrifugal 
term corresponding to an "angular momentum" m = ±1. 

One can develop an S matrix formalism in this representation, and extract the 
L-Z rotation from the branching properties of the solution at = 0: indeed at the 
level crossing the frequency of the motion inverts its sign. It should be further pointed 
out that, as long as Z is quadratic in the spinor amplitudes, its branching behavior is 
accordingly modified with respect to the spinor case. In Appendix II we give the main 
steps of the solution. 

Before going to numerical results, let us summarize the qualitative features of the 
motion: 

a) in the "normal" regime the representative undergoes uniaxial rotations (precession 
with time-dependent frequency having an average value of ±2i?o/7r); b) in the L-Z 
regime a discrete rotation, involving the population as well as the coherences, occurs; 
c) in the next "normal" regime the bias ci;o(^) has opposite sign with respect to regime 
(a), so that the sign of precession is inverted in the L-Z region. 
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Let us consider the following 3 situations: T fa 1, 0, —1. 

In the first case we have {y 1), so that, after crossing, an initially "magnetized" 
state is preserved, apart from a small correction; in the second case, where the transition 
probability equals 1/2, it is turned into a "fully unmagnetized" state; in the third case 
(i/ 3> 1), magnetization reversal occurs. 

Fig. 1 is the three-dimensional plot of the representative v. The interlevel couphng 
A is of the order of the frequency Qq- 

The L-Z crossings can be found in the regions where the trajectory inverts its path. 
Furthermore, at the crossings the vector v undergoes a discrete rotation {9, 0}. The 
vertical angle 9 is determined by the transition probability between the two levels, which 
goes as exp{—2n ■ u}: here, where we started with Z — —1, one can see that \Z\ is 
reduced accordingly at each crossing. 

In Fig. 2, Z{t) is plotted over 10 cycles of the external field. The fast oscillations 
correspond to the "normal" regime, where Z is on average constant; the "normal" 
regimes are separated by kinks, produced by the crossings. 

Here we have 9 ^ 7r/3, as three kinks are needed to reach the complete inversion of 
Z. The population Z (as well as the coherences X, Y) undergoes large scale oscillations, 
following the bias frequency. A planar {X, Y) plot would exhibit similar oscillations 
associated with (p. 

One could argue that a zero- average oscillating bias Bz = uJo{t) must generate 
a zero-average Z(t). On the contrary, "single spin magnetization" is possible, with 
hamiltonian dynamics. Indeed, by operating on the "normal" evolution through Qq, we 
obtained trajectories where Z{t) keeps a fixed sign. 

This can be understood as follows. Although the L-Z rotations always act additively 
on V, if the interposed "normal" evolution inverts the planar components, the sequence 
of the 9 rotations acts on Z{t) alternatingly rather than additively. 

Under these conditions if, e.g., one starts with Z{0) — —1, and 9 = tt/S, at the 
first crossing Z is raised at Z 0.3, but at the next crossing, it turns back to Z — — 1. 
An example is shown in Figs. 3, 4, 5, referring to Z{t), v(t), Z{B) respectively. 

In this non-adiabatic regime in general we found strong irreversibility; the 
"symmetry-breaking" solution instead is almost perfectly reversible. This is made clear 
in Fig. 5, where practically no hysteretic effects arc found. 

When the transition probability equals |, the population, when started at Z — —1, 
turns into a perfectly balanced superposition of the two levels {Z — 0) and the 
vertical angle is: 9 — 7r/2. The "in-phase" evolution follows the sequence Z — 
— 1 ; ; 1 ; ; —1..., the "out-of-phase" evolution is given hy Z — — 1 ; ; — 1 ; 0.... The 
latter case again generates a completely reversible process. 

At much larger values of the transition probability the magnetization aligns to the 
bias at each step. The field evolves slowly with respect to the two-level system: we are 
in the adiabatic regime, and we obtain very small deviations from reversibility in the 
plots Z — Z{B) (not shown here). 



The Landau-Zener Model with Decoherence: The case S=l/2 8 

5. Irreversible dynamics: homogeneous case 



The simplest non-hermitean generalization of the previous model includes population 
damping and decoherence. This leads to an homogeneous equation for v: 



X + 7i X = -ujqY 

Y + -f^Y =uoX - AZ (11) 
Z + 73Z = AF 

In compact form, we have = ■ v. Notice that the isotropic case 71 = 72 = 73, 
apart from a Poincare sphere with decaying radius, is identical to the unitary case. 

We thus consider the general case. The large time asymptotics of the solution of 
Eq.|TT| can be determined from the eigenvalues Pjif) of the operator B{t) [jl5[. 

We limit ourselves to few preliminary results, for the sake of simplicity in the 
uniaxial case: 71 = 72 = 7r; 73 = 7- 

The equation for the right eigenvectors Ej(t) reads: 

5(t)-E,(t) = p,(t)-E,(t) (12) 

It can be shown that the large time asymptotics is determined by the following 
fundamental solutions: 

\j{t) ^ exp{i ■ I Pj{u) ■ du] ■ Ej(t). (13) 

Since here we have j = 1, 2, 3, the eigenvalues can be analytically computed. 
It is convenient to shift the variable p : p y = p + •jr + [5 = (7 — 7r)/3]. The 
eigenvalues are then: 

Hi = u + V 



1/2,3 = 2 ' 



u-v)] (14) 

where are the roots of the equation: 

t'+ft-g^ = (15) 



f = 26[m' + 6' - ^A'] 
g = ^[{Qtr + A' - 36'] 



In the range of parameters we are interested in, both roots of eq. are real- valued. 
We keep terms up to the order 0(l/t^) in the expansion in powers of t, and get: 

1 

Pi{t) ~ - 7 -27c 



(nt) 



2 



P2,3(t)~ -7r + 7c^±^[^^t + y] (16) 

7c = S{S' - '-A') 
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The first eigenvalue gives the correction to the damping 7, while ^2,3 in their real 
part, give the correction to 7,.. It is worth noticing that these corrections act on 7 
and7r with opposite signs: in other words, decoherence is slowed down at the expense 
of damping, and viceversa, according with the sign of 7c. 

The imaginary parts of P2,3(^), related with the oscillating behavior and with the 
branching properties around t = 0, are not influenced by decoherence and damping. 
One verifies that, apart from a factor of 2, which is expected in going from the spinor 
amplitude to the vector amplitude, the result coincides with the semiclassical momentum 
pit), given after eq. [l^. 

In the example of Fig. 6 we assume pure decoherence in the adiabatic regime: the 
population Z follows the bias. Notice that indeed Z is damped as expected from the 
correction to 7 computed previously(see eq. |16|, with 7c > 0). 

The values of X and F, between crossings, sustain an oscillating Z with an {Z) 
sensitively different from zero. 

The symmetry breaking solution, found in the previous section, survives to 
decoherence, as shown in Fig. 7. Here we have the conditions of Fig. 3, but the 
decoherence parameters of Fig. 6. One can realize a "magnetized" single spin, by 
means of an oscillating bias. Here the magnetization is in the direction of the bias 
field; the situation discussed in ref. is rather different: there the magnetization is 
orthogonal to the bias and Z relaxes toward the ground state. 

When 7c < 0, in the adiabatic regime, the dominant effect is that the population 
decay is enhanced and Z{t) rapidly evolves towards zero. 



6. Inhomogeneous case: relaxation from upper to lower level 

The incoherent processes in a two-level system include internal transitions (i.e. interlevel 
transitions mediated by the surroundings). 

Here we analize the effects of incoherent relaxation from the upper to the lower 
level. 

At the L-Z crossings the ground state flips between the two levels, and the 
inhomogeneous term must change its sign accordingly; we are then led to the 
inhomogeneous equation: 



X + 71 X = -ujQ Y 

y +72F = cuoX - AZ (17) 
Z +73(Z + ^) = Ay 

If one were to choose a sign for the inhomogeneous term in Eq. [l^, one would 
describe a situation in which the system would be forced to decay towards a given level. 
This is appropriate only provided that ujq has a fixed sign. 
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The calculations presented in Fig. 8, 9 refer to the adiabatic, strongly damped 
regime. The population Z starts switching from a "magnetized" state to the opposite 
one, exactly when the external field changes its sign. 

The escape from the instability is markedly different from the relaxation towards 
the stable ground state (see Fig. 9), but this is not the origin of the hysteretic behavior 
shown in Fig. 8. 

In the approach towards a fully "magnetized" state, the coupling between Z and 
the coherences X, Y plays a relevant role (see the oscillanting branches in the hysteresis 
cycle). When the bias retraces its path backwards, the system is already in a fully 
"magnetized" state Z = ±1, X = Y = 0. Only the onset of the instability, as the bias 
changes its sign, is able to put the system in motion again. 

In conclusion, we have two regimes (incoherent X, Y ^ and coherent X, Y ^ 
respectively) both compatible with the same value of the bias. 

So far we discussed a situation where the period 27t/Qq is much larger than the 
decay times of the two-level system, so that the external field can be considered as 
quasi-static. As the decay time gets longer, the system does not relax fast enough 
towards the "fully magnetized" state: an hysteretic behavior survives, as shown in Fig. 
10. 

We conclude this section with a last example: we added the relaxation towards the 
ground state to the system discussed in Fig. 7; the frequency Qq is comparable with A. 

In Fig. 11, the lastly oscillating portions refer to the "normal" evolution, and the 
kinks to the crossings. Notice that the oscillating parts are attracted towards the "fully 
magnetized" state, which is never reached. The system shows an hysteretic limit cycle 
(see Fig. 12). 



7. Conclusions 



We examined the time evolution of a two-level system in the presence of an oscillating 
external field, both in the unitary case and with decoherence and dissipation. We 
solved, in this context, the Landau-Zener problem in the density matrix formalism, and 
showed that, in the frequency representation, the population Zi^uj) satisfies a Schrodinger 
equation for an inverted harmonic oscillator in D=2. 

We indicated a procedure allowing to interpret the various dynamical regimes by 
means of few intuitive rules. We found that it is possible to sustain a magnetization 
parallel to the bias for the single spin. 

We became recently aware of related work |0, where instead the spin is magnetized 
in the direction orthogonal to the bias; furthermore, while our result holds in the 
homogeneous case, in Ref. relaxation towards a single level is assumed. 

We have also shown that one can explicitly solve the Landau-Zener problem in the 
non-unitary case, by means of a suitable extension of WKB methods. We stress that 
in our approach the parameters of the equation are completely arbitrary, so that no 
"adiabatic elimination" of fast variables is needed ||14| . 
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We find that the corrections to 7 and 7^ have opposite sign. The branching 
properties of the solution are not modified with respect to the unitary case: the 
scattering matrix is again a function of the single L-Z parameter v. 

When decay from the upper to the lower level is taken into account, we find 
hysteretic behavior: this effect originates in the coexistence of a fully magnetized state, 
where the coherences X, Y are strictly zero, and a partly magnetized state, with nonzero 
coherences. This interpretation holds in the overdamped case; at slower dampings an 
hysteretic limit cycle survives, related with non-adiabatic effects. 

8. Appendix 1. Inequalities following from complete positivity 

The conditions for complete positivity of the evolution, when written for matrix A^, 
are: 

< 71 < 72 + 73 < 72 < 73 + 71 < 73 < 71 + 72 

together with: 

^{l' + \cl) < 7? -(72 -73)^ 
^{(^'+\CI) < 72^-(73-7if 

4/3(a7 - \c^Cr) - 2C2{\aCr + ]^iC^) > (71 + 73 - 12){P' + l^,) + 



(72 + 7i - 73) [tt^ + 4 ~ 4 ^3 - ^ (72 - liY] + (73 + 72 - 7i)(7^ + ^ Ci) 

In our problem we further must take into account a;o(^), which defines the scale of times, 
and A. Thus A is the single relevant parameter of the model. 



9. Appendix II. Solution of the Landau-Zener problem from the master 
equation 

In this Appendix we summarize the main points of the solution of the L-Z problem, as 
obtained by starting from the master equation; we omit those details that can be easily 
implemented from the standard solution. 

In the non-dissipative case the Landau-Zener problem is written in the form: 



The Landau-Zener Model with Decoherence: 



The case S=l/2 



12 



x = - iytt) ■ Y 

Y = (nt) -X -A- Z (18) 
Z = A-Y 

The variable Y can be eliminated from the third equation, and one obtains: 

Z + A"^ ■ Z = Ant ■ X (19) 

n 

x= - —t-z 

A 

Let us define the new variable (p{t) = X(t) + £t • Z(t); it is readily verified, by 
inspecting the second equation in system |19|, that from (j){t) one can recover the entire 
solution: Z{t) = (A/n) ■ ^(t); X{t) = <p{t) - 1 ■ <p{t) . 

The Fourier transform 0(c<j) = exp{—iujt} ■ (j){t) ■ dt solves the equation: 

H4){uj) = 

,A^ . d 



(20) 



Let us write H as: 



H = h-uo 

We find that the function x{^) = ^ ' 0(^)) which, apart from a constant factor, 
is the Fourier transform of the population Z{t) : (xi^) = £ ' ^i^))^ satisfies the 



Whittaker equation 18 



hxiuj) = 0. (22) 



Notice that eq. ^ has the form of a Schrodinger equation for the inverted harmonic 
oscillator in D=2, with "angular momentum" m = ±1; (clearly in our case u varies 
over the whole axis). 

The equation can be written in the canonical form of the hypergeometric confluent 
equation (HCE) by letting: q = cu'^ ; x{^) = ^^P{'^} ' {qY^"^ ' u{q). 

In terms of the new variable q '■ q = (i^) ■ ^, one finally gets, for the function 
w{0 = u{q): 

^ ■ + {2-i)-^^w-{l + iv)-w = Q (23) 
_ A2 

Two linearly independent solutions, respectively regular at infinity and around zero. 



are: 



,2 n, ,2 



X\iy>) = expj— — } ■ ■ \E'(a, 2 
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X2H = exp{— } • a; ■ $(a, 2, 



a = 1 + ii/ 

where $, ^ arc the hypergeometric confluent functions of flrst and second type. 

In the L-Z problem the large time region is characterized by high frequencies, since 
there the level spacing always dominates over the interlevel coupling. Hence the solution 
Xi(a;), having a well-defined behavior at infinity, is the natural choice in our case. 

We consider the behavior of the function ^(a, c, z): 

*(a,c,^)«(^)-«-[l + 0(-)], (^»1), 

z 

and recall that ^ satisfies the identity: 

-^{a,c,z) = (2;)^"'' • *(a - c + 1,2 - c,2;). 

If one takes the former properties into account, and uses the identity: (t){oj) — 
a;)), which follows from being real, one obtains: 

+ ^-exK^}. ^. 0.—) 

(24) 

One recovers the behavior of the population from Z{ijj) = i^u) ■ (f){(jj); we have, as 
a; — >• 00, Z{u!) f» l/u. 

More precisely, apart from constant factors: 

Z{u;) = z^-exp{-nu/2}iexp{z^}i^)-'-+exp{-z^}i^^ 

It is easily verified that the correct transition probabilities can be extracted from 
the branching properties around a; = of the solution Z{u!). 

We remind that the standard treatment refers to spinor amplitudes, while Z{u!) is 
quadratic in such amplitudes: this explains the factor of 2 in the exponent of ou. 

The coherence X[u}), can be found from the identity: 

X(a;) = 0H + ^M(a;)). 

We add the large u; behavior of X{u;): 

X{u;)^^{l-2{^r)exp{-7ri.)/2}[exp{t^}i^^ 

In summary, while the Schrodinger equation leads to the inverted harmonic 
oscillator in D=l (Weber equation) in the representation of times, the master equation 
leads, for Z{uj), to the D=2 inverted harmonic oscillator. Hence, in the representation 
of frequencies, it is possible to formulate and solve the L-Z problem in terms of the S 
matrix. 
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Figure 1. A typical trajectory v(t),v = X,Y,Z, starting at = 1. The "normal" 
evolution corresponds to the fast uniaxial rotations: here four branches of this sort can 
be identified. At the crossings the trajectory inverts its path and makes a transition 
to a new branch. We have A = 0.01, Oq = 0.02; here and in all other figures Bq = 1. 
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Figure 2. Time evolution of Z over 10 cycles of the bias field. Three kinks are 
needed in order to have a complete inversion of the magnetization: hence 6 = ir/3. 
Here A = 0.12, Qq = 0.063. Here and in all the Z(t) plots the scale of time is in 
arbitrary units. 
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Figure 3. Time evolution of Z over 10 cycles of the bias field, with A = 0.12, 
Oo = 0.0682. With this choice, at each crossing the vertical rotation changes its sign: 
the result is an explicit symmetry breaking: < Z{t) >< 0. 
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Figure 4. Trajectory v(t), parameters as in Fig. 3. 
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Figure 5. Plot oi Z = Z{B), from the solution of Fig. 3. Notice that, in spite 
of being far from the adiabatic regime, the system shows very slight deviations from 
complete reversibility. These deviations are very sensitive when the system is tuned 
out of the particular regime displayed in Figs. 3, 4, 5. 
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Figure 6. Plot of Z[B) for a single cycle of the bias, with initial condition 
Z{t = 0) = l.The parameters are A = 0.3, = 0.033, 7^ = 0.01,7 = 0. 
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Figure 7. Survival of magnetization in the presence of decoherence: plot of Z{t) 
over ten periods of the bias, with A = 0.12, Oq = 0.0682, 7^ = 0.01,7 = 
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Figure 8. Approach to the hysteretic Umit cycle in the strongly damped case: Z{B) 
with -fr = 0.035, 7 = 0.07 and A = 0.05, = 0.02. Notice that Z is almost 
perfectly constant in the fully magnetized branches, where the coherences X, Y are 
zero. The magnetization reversal starts at the zeros of the bias field. 
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Figure 9. Plot of the magnetization Z{t) in the same conditions of Fig. 6a: the cusps 
corresponding to the onsets of magnetization reversal separate perfectly identical, but 
inverted, kinks. The result is a sequence of pulses lacking specular symmetry. 
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Figure 10. Same as in Fig. 6a, but with smaller damping: 7^ = 0.01, 7 = 0.02. 
Here the complete magnetization is never reached, because the relaxation towards the 
ground state has a time scale comparable with the period of the bias. An hysteresis 
is still obtained; the behavior of Z{t), not displayed here, is again a sequence of 
asymmetric pulses. 
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Figure 11. Plot of Z{t) in the non-adiabatic regime; the parameters are as in Fig. 5, 
but now we added the relaxation towards ground state with 7 = 0.02. It is instructive 

to compare this result with the hamiltonian case of Fig. 3b, where we had the same 
values of A and ^q: there Z never changes its sign, here it is attracted towards the 
lower level. Again each pulse is strongly asymmetric. 




Figure 12. Approach to the hysteretic hmit cycle of the solution discussed in Fig. 8a: 
one can identify the cusps at the onset of magnetization reversal. 



